rm(list = ls())

data <- read.csv("govdata.csv")


gov.nodecs <- lm(inc.twpct ~ presgov.twpct + govvote.lag +  propertypercapl.6mo +   medincK + as.factor(countyfips) + as.factor(year), data = data)

gov.decs <- lm(inc.twpct ~ presgov.twpct + govvote.lag +  propertypercapl.6mo +  disdecs.all.6mo + medincK + as.factor(countyfips) + as.factor(year), data = data)

gov.tds <- lm(inc.twpct ~ presgov.twpct + govvote.lag +  propertypercapl.6mo + disdecs.all.6mo + turndowns.6mo + medincK + as.factor(countyfips) + as.factor(year), data = data)